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ACTIVE  AND  PASSIVE  REMOTE  SENSING  OF  ICE 


Principal  Investigator:  Jin  Au  Kong 

Sell -Annual  Progress  Report 

This  Is  a  report  on  the  progress  that  has  been  made  In  the  study  of 
active  and  passive  remote  sensing  of  Ice  under  the  sponsorship  of  ONR  Contract 
N00014-83-K-0258  during  the  period  of  February  16,  1983  -  July  31,  1983. 

During  this  period  we  have:  (1)  derived  the  dyadic  Green's  function  for  a 
two-layer  anisotropic  medium;  (2)  derived  the  backscattering  cross  sections 
and  the  bistatic  scattering  coefficients  for  a  two-layer  anisotropic  random 
medium;  and  (3)  participated  in  a  planning  meeting  at  the  Cold  Regions 
Research  and  Engineering  Laboratory  (CRREL)  concerning  the  winter  microwave 
remote  sensing  measurements. 

The  dyadic  Green's  function  for  a  two-layer  anisotropic  medium  has  been 
derived.  The  anisotropic  medium  Is  assumed  to  be  tilted  uniaxial.  Inside  the 
uniaxial  medium  there  are  two  characteristic  waves,  l.e.  an  ordinary -wave  and 
an  extraordinary-wave.  An  Incident  wave  with  either  TE  or  TM  polarization 
produces  both  TE  and  TM  waves  when  It  Is  reflected  at  the  Interface  of  an 
Isotropic  and  a  uniaxial  anisotropic  media  because  of  Its  coupling  Into  both 
the  o-wave  and  the  e-wave.  With  the  availability  of  the  dyadic  Green's 
function,  mapy  electromagnetic  wave  radiation  and  scattering  problems  In 
layered  uniaxial  media  can  be  solved.  A  manuscript  has  been  prepared  for 
submission  to  a  journal  for  publication  [Appendix]. 
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The  backscatterlng  cross  sections  and  the  bl static  scattering 
coefficients  for  a  two-layer  anisotropic  random  medium  have  been  derived. 

The  Bom  approximation  Is  used  along  with  the  dyadic  Green's  function  for  the 
two-layer  anisotropic  medium  to  calculate  the  scattered  fields.  The 
theoretical  results  are  used  to  match  published  experimental  data.  We  are  now 
for  the  first  time  able  to  compute  contributions  from  both  the  background 
anisotropy  and  the  anisotropy  of  the  permittivity  fluctuations.  A  manuscript 
Is  being  prepared  to  document  the  theory  and  the  theoretical  calculations. 

We  have  participated  at  the  planning  meeting  at  CRREL  concerning  the 
winter  microwave  remote  sensing  measurements.  On  June  17,  1983,  the  Principal 
Investigator  went  to  a  meeting  at  CRREL  concerning  the  measurement  of  the 
physical  and  electrical  properties  of  fresh  water  Ice  and  artificially  grown 
sea  Ice  for  the  determination  of  the  effects  of  these  properties  on  the 
emlsslvlty  and  reflectivity  of  the  Ice.  A  presentation  was  made  to  stress  the 
Importance  of  measuring  the  cross-sectional  profile  of  the  Ice  which  has  been 
shown  to  play  a  very  significant  role  In  the  random  medium  model. 
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APPENDIX 


<ad1c  Green's  Functions  for  Layered  Anisotropic  Medium 


J.  K.  Lee  and  J.  A.  Kong 

Department  of  Electrical  Engineering  and 
Computer  Science  and  Research  Laboratory 
of  Electronics 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


Abstract: 


The  dyadic  Green's  functions*  (DGF)  for  unbounded  and  layered  anisotropic 
media  have  been  obtained.  The  anisotropic  medium  is  assumed  to  be  tilted 
uniaxial.  With  the  availability  of  the  DGF's,  many  problems  involving  radia¬ 
tion  and  scattering  of  electromagnetic  waves  can  readily  be  solved. 


N  *-  A 


I.  Introduction 


The  dyadic  Green's  functions  (DGF)  for  various  configurations  have  been 
studied  extensively  in  the  literature^"^.  In  applications  to  microwave 
remote  sensing,  the  problems  of  radiation  and  scattering  of  electromagnetic 
waves  in  layered  uniaxial  media  have  presented  a  fundamental  difficulty  in 
the  absence  of  a  useful  dyadic  Green's  function.  The  DGF  for  a  layered  uni¬ 
axial  medium  with  the  optic  axes  perpendicular  to  the  planes  of  stratification 
has  been  derived^.  However,  in  most  practical  applications,  the  optic  axes 
of  the  uniaxial  medium  is  tilted.  It  is  the  purpose  of  this  paper  to  present 
the  results  of  the  DGF  for  both  an  unbounded  and  a  layered  anisotropic  medium 
corresponding  to  a  uniaxial  medium  with  a  tilted  optic  axis.  In  order  to  cast 
the  final  results. in  a  simple  and  i nterpretabl e  form,  the  reflection  and  trans¬ 
mission  coefficients  at  the  boundary  surface  of  the  anisotropic  medium  are  also 
obtained. 
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II.  Dyadic  green's  Function 


In  order  to  facilitate  the  derivation  of  the  dyadic  Green's  function  (DGF) 
for  an  anisotropic-medium  layer,  we  first  consider  an  unbounded  anisotropic 
medium  whose  permittivity  tensor,  c,  is  taken  to  be  uniaxial  with  its  optic 
axis  tilted  off  the  z-axis  by  an  angle  y  in  the  y-z  plane  as  shown  in 
Fig.  1.  The  permittivity  tensors  and  e,  before  and  after  the  tilting, 

respectively,  are  as  follows: 


fcll 

S 

e  *  0 


e22  e23 


32  e33 


where 


«n  *  e. 


622  *  e  cos2’*'  +  ez  s‘*n2,*'» 


e23  *  e32  "  ^£z  "  cos  *  s^n 


e33  ■  e  s1n2ij>  +  &z  cos2i 


■  ,» 1 
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The  DGF  for  this  anisotropic  medium,  §(r,  r* ) ,  satisfies  the  vector  wave 
equation 

v  x  v  x  §(?,  r* )  -  u)2pe  •  B(r,  r ‘ )  =  !  6 ( r  -  r ' )  (4) 

where  w  is  the  angular  frequency  and  y  is  the  isotropic  permeability  of 
the  medium.  Me  write  the  Fourier  transform  pair  of  DGF  as  follows 

S(r,  ?')  *  — - —  [  d3E  G(fc,  r* )  e1^  *  ^  (5) 

(2ir)3  ■*  -«° 

S(E.  f' )  «  f "  d3r  8(r,  ?' )  e'1f  '  (6) 

}  .00 

Substituting  (5)  into  (4)  and  using  the  identity 

6(r  -  r1)  *  f"  d3R  eiR  *  (7) 

{  2ir )  3 

we  obtain 

G(ic,  ?')  *  -[E  •  k  +  u^ye]*1  e'lR  ’  r‘  (8) 


where 

ro 


(9) 


Performing  the  matrix  inversion,  we  obtain 


S(«c,  ?•) 


where 


,e-iE  *  ?' 

«2u  D^(R)  02(R) 


*11 

A12 

A13 

^1 

^22 

^23 

A32 

A33 

(10) 


V*) 

»  k  2 
*x 

+  V 

+  kz 

,2  -  (D2ye 

(11a) 

02(C) 

*  £kx 

1  +  e22 

k  2 

y 

+  £33kz2 

+  e232kyk; 

Z  "  “2»*ee2 

(lib) 

A11  * 

kx2°l 

+  k/“ 

2u(e 

*  £z>  - 

«2uD2 

(12a) 

H 

CM 

« T 

*21  * 

kxkyDl 

+  01 

2ykx(e  - 

ez)  cos  V 

(k  cos  <p  - 

kz  sin  \p) 

(12b) 

*13  * 

*31  * 

kxkzDl 

+  (1) 

2ukxuz  - 

e)  sin  tp 

(k  cos  p  - 

kz  sin  ij>) 

(12c) 

*22  * 

k/D, 

-  <D2uCekx2 

+  e33ky2 

+  e33kz2 

’  “3yee33^ 

(12d) 

*23* 

*32  * 

kykz°l 

+  ID 

2w23(ky2 

+  kz2  -  a 

>2ye) 

(12e) 

*33  * 

kz2Dl 

-  u2y[skx2 

+  e22kyZ 

+  e  k  2 
e22Kz 

-  a)2uee22] 

(12f) 

In  the  derivation  of  (11)-(12)  we  made  use  of  the  relationships 


e22e33  "  e23  =  e£z  ^13a^ 

e22  +  e33  =  e  +  V  (13b) 

Substituting  (11)-(12)  into  (5),  we  perform  the  integration  over  kz>  The 
poles  of  the  integrand  occur  at  the  zeros  of  (R)  and  D2(£)  denoted  by 
±kz°  and  kz  3  k*u  and  k®d,  where 


,°  *  0*»2ve  -  kx2  -  ky2]1/2 

(14a) 

eu 

z 

ed  j 

k  ’  "  23  ky  1  1  ^e23ky2  ”  elle33kx2 

"  E22e33ky2  +  “2^eze33]1/2. 

*  J 

e33  e33 

(14b) 

Here,  kz°  (or  -kz°)  and  k®u  (or  k®d)  are  the  z-components  of  the  up¬ 
ward  (or  downward)  propagation  vectors  for  an  ordinary  wave  and  an  extraordinary 
wave,  respectively.  We  also  notice  that  the  dispersion  relations,  D-j ( R )  =  0 
and  D2(R)  3  0,  can  be  obtained  by  simply  transforming  the  R-variables  in  the 
dispersion  relations  of  vertically  uniaxial  medium^,  i.e.  by  replacing  k  , 

A 

k  ,  ^  k  ,  k  cos  -  k  sin  ip ,  k  sin  ip  +  k  cos  \p  in  (All)  and  (A12) 

Jr  i  x  y  £  y  £ 

of  Ref.  [2]  which  are  shown  here: 


(18c) 


B22  3  kx2  cos2,<'  (18d) 
B23  *  B32  *  'kx2  sin  *  cos  *  (18e) 
B33  s  kx2  sin2^  (18f) 
C]1  =  kx2(ky  sin  ip  +  k®u  cos  ip)2  (19a) 


C-|2  *  Cg-j  =  kx(ky  sin  +  k2  cos  <p)[ky(ky  sin  ip  +  kz  cos  -  u)2ue  sin  t|>] 

09b) 

C13  3  C31  *  kx^ky  S1n  *  +  kiz  cos  *)[k®U(k  sin  ip  +  k®u  cos  ip)  -  u>2ue  cos  <|>] 

(19c) 

C22  *  Cky(ky  sin  *  +  kzU  cos  ^  “  u2ue  sin  (19d) 

C23  3  ^32  =  Cky(ky  sin  *  +  kzU  cos  ^  ”  w2ye  sin  ^C^^y  sin  ♦  +  kzU  cos  ^ 

-  w2ue  cos  ip]  (I9e) 

C33  *  [k®u(ky  sin  ip  +  k®u  cos  ip)  -  w2ue  cos  ip]2.  (19f) 

The  following  relations  are  used  in  the  derivation  of  (16)-(19). 

e[k  2  +  (k  cos  ip  -  k  e  sin  ip)2]  *  -e  [(k  sin  *  +  k  6  cos  \ p)2  -  o>2ue]  (20a) 

X  y  2  z  y  4 

eCRx2  +  ky2  +  k®2  -  ai >2uez]  =  (e  -  ez)(ky  sin  ip  +  kze  cos  ip)2 


(20b) 


(21a) 


D2(kz  -  k2  )  *  (e  -  ez)[kx2  +  (ky  cos  <p  -  kz  sin  <p)2] 


Di(kz  =  k  e)  * - -  [(k  sin  tp  +  kze  cos  ip)2  -  w2ue]  (21b) 


e  , eu  i ed 

:  =  kz  or  kz 

22  si*1  ♦  +  e23  cos  ^  s  ez  sin  ^ 

(22a) 

: 2 ^  sin  ip  +  cos  ip  =  ez  cos  ip 

(22b) 

22  cos  tp  -  e23  sin  ip  =  e  cos  ip 

(22c) 

22  cos  ip  -  £^2  cos  ip  *  -e  sin  ip 

(22d) 

;  -  e22  =  (e  -  ez)  si n2ip 

(23a) 

2  -  6 22  =  (ez  "  e)  cos2^ 

(23b) 

-  £33  =  (e  -  £z)  COS2ip 

(23c) 

z  ■  E33  *  <Ez  •  e)  s1n2'" 

(23d) 

For  z  <  z',  we  obtain  !(?,  r* )  in  a  similar  manner  by  noticing  that 
Im(-kz°)  <  0  and  Im(k®d)  <  0. 

Finally  we  cast  the  dyadic  Green's  function  for  the  unbounded  anisotropic 
medium,  S(r,  ?' )  in  the  following  form: 
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<>  •  ftIt  (25f) 


R°  •  E  -  k,°i 

P  z 


R •  -  R  +  kfdz 
p  z 

E  •  e  •  Ee 


Ku  -  e  •  Ke 


Cp  *  k„X  ♦  y 

A  A  A 

z'  *  y  sin  ip  +  z  cos 


We  notice  that  o(±kz°)  is  a  unit  vector  in  the  direction  of  the  electric 
field  for  an  ordinary  wave  and  e(k®u)  or  e(k®d)  is  a  unit  vector  in  the 
direction  of  the  electric  field  for  an  extraordinary  wave.  It  is  seen  that 

a 

o  is  linearly  polarized  perpendicular  to  the  plane  formed  by  the  optic  axis 
(z'-axis)  and  the  propagation  vector  (E°-vector)  as  it  represents  an  o-wave, 

A 

and  e  is  linearly  polarized  parallel  to  the  plane  formed  by  z'-axis  and 
Re-vector  as  it  is  an  e-wave^.  We  also  notice  that  for  an  o-wave,  the  di¬ 
rection  of  power  flow  determined  by  the  Poynting  vector  E  x  R*  is  the  same 
as  the  direction  of  propagation  vector  (R°),  while  for  an  e-wave  the  direc¬ 
tion  of  power  flow  is  different  from  that  of  propagation  vector  (Re).  In 
fact,  for  an  upward  propagating  e-wave,  the  electric  field  (E),  the  mag¬ 
netic  field  (fi)  and  the  Poynting  vector  (E  *  H*)  are  in  the  direction  of 


e(kze) ,  o(kze)  and  Eu>  respectively. 

The  result  checks  with  the  expression  of  DGF  for  a  vertically  uniaxial 
medium^  in  the  limiting  case  as  y  »  0.  For  an  isotropic  medium  where  <i>  *  0 
and  e  31  e2,  the  o-wave  becomes  a  horizontally  polarized  TE  wave  and  the  e- 
wave  becomes  a  vertically  polarized  TM  wave  as  follows. 


o(±k  °)  -*■  h(±k  )  *  —  z  x  E 

2  z  k 

(26a) 

p 

e(k*u)  -  v(k2)  =  1  h(k2).  x  E 

(26b) 

e(k®d)  -  v(-kz)  *  1  h(-k2)  x  R 

(26c) 

where 


kp  ’  ("x2  +  V>,/2 

(27a) 

kz  >  <k*  -  koM1/2 

(27b) 

k  a'  un/yeT 

(27c) 

E  *  E  +  k  z 

P  2 

(27d) 

R  *  E  -  k  z. 

P  2 

(27e) 

(27e) 
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Next  we  consider  the  two-layered  stratified  medium  as  shown  in  Fig.  2. 
The  upper  and  lower  media  are  isotropic  and  characterized  by  eq  and  e2> 
respectively.  The  medium  in  the  middle  is  anisotropic  as  described  in  the 
previous  section.  In  the  principal  coordinate  system  it  is  characterized  by 
ej  and  e^z-  The  permittivity  tensor  after  the  tilting  of  the  optic  axis 
is  given  by 


where 


ell  *  el *  (29a) 
e22  *  el  C0s2,<'  +  elz  sin2,<'  (29b) 
e23  “  e32  *  (elz  "  £i)  sin  cos  *  (29c) 
e33  "  ei  sin2^  +  Eiz  cos2'!'-  (29d) 


All  three  regions  are  assumed  to  have  the  same  permeability  y.  Let  the  source 
be  in  region  0.  DGF's  S0g(r,  ?'),  l10(r,  and  S20(r,  r’)  satisfy 


V  x  7  x  ^(r,  ?')  -  o»2ye0  ^{r,  ?')  *  !  5(r  -  ?'),  z  >.  0  (30a) 


7  x  7  X  r,  r')  -  w2ye1  •  !L,0(r,  r')  3  0, 
7  X  7  X  S^r,  )  "  «2W£2  ^20^*  =  0» 


-d  <  z  <  0  (30b) 


z  <  -d  (30c) 


and  the  appropriate  boundary  conditions  at  z  =  0  and  z  *  -d  and  radiation 
conditions  at  z  *  ±«.  The  first  subscript  of  the  DGF  refers  to  the  region  of 
the  field  and  the  second  subscript  to  the  region  containing  the  source. 

Since  the  DGF  for  an  Infinite  unbounded  media  assumes  the  form  in  (24), 
it  is  appropriate  to  assume  the  DGF's  for  this  two-layered  geometry  in  the 


following  form: 


®oo<f* ?,)  ■rr[L*Ar-'[^(-koz>  •  0  r  +  "hh  iio(koZ)  e  0 


iR_  •  r 


+  "HV  ;0(lt02)  *  "  r]  *iO(-kOZ)  +  [;o(-k02>  *' 


lie  •  r  «  ilc  •  r  A  -1R  •  r1 

♦  “vv  vo<koz>  e  *  RVH  h(koz>  e  3  *0(-K0Z))  e 


0  <  z  <  z' 


«  iR  0  •  r 

«!0<?-  ^  *  “T  f  dkxdky  ~  e  1  +  Bun  o(k?J  e 


(31a) 


lie,0  •  r 
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i 

! 


¥ 


Cil 


■  »H  -  eu  lE,e*r  . 

*  *He  e  +  BHe  e*klj*  e  ^  ho*'koz' 


II!,0.  f 


*  [\o  °<-kU>  e  +  BVo  o(klz>  e 


.  IE,0-? 

o  \  „  1 


.  iR,®  *  r 


♦  \.  °<kn> 0 


:„.eu,  JEie-'r 


+  BVe  e(k=“)  e 


-iR  •  r* 

]  vo<-koz>  0  • 


-d  <  z  <  0 


(31b) 


Wr*  r*  ^  *  a!2  1L.  J  M_k2z)  +  *HV  v2^’k2z^  V‘koz) 

°Z  1 

.  .  .  iR,  •  r  -IK  •  r‘ 

+  C*vv  v2("k2z^  +  XVH  h2^“k2z^  vo^“kozM  e  e  * 


z  <  -d 


(31c) 


where 


Mtkiz>  *7*  “  ki> 

P 

*i(kiz>  •  r- *’i(kiz)  x  lr 
Ki 


v1(-k1z)  *  -h1(-k1zJ  x  Ri' 
*1 

-  o  z'  *  *1° 

°(kiz)  s  7: — ttt 

12  | z 1  x  E1 °| 


1  =  0,  2 


i  *  0,  2 


i  *  0,  2 


(32a) 


(32b) 


(32c) 


°Kz> 


z1  x 

li*  ><  *i°l 


17 


kj  ■  w/ye^ ,  i  =  0,  1,2. 


The  boundary  conditions  to  be  satisfied  at  z  *  0  and  z  *  -d  are 


i  «  ^(f.  *') 

*  z  x  ^(r,  r') , 

z  =  0 

(33a) 

z  x  G10(r,  r') 

*  z  x  (^(r,  r* ) , 

z  -  -d 

(33b) 

Z  x  7  x  Gqq(?, 

*  * x  7  * 

z  -  0 

(33c) 

Z  x  7  x  §.jq(?. 

f)  -  Z  x  7  x  ^(f,  ?) 

z  =  -d 

(33d) 

where  (33a)  and  (33b)  refer  to  the  continuity  of  tangential  electric  fields 
and  (33c)  and  (33d)  refer  to  the  continuity  of  tangential  magnetic  fields  at 
each  Interface.  Substituting  (31)  into  (33),  we  obtain  eight  linear  algebraic 
equations  for  the  eight  unknown  coefficients.  However,  it  is  very  difficult 
to  solve  these  equations  analytically  and  to  express  the  unknowns  in  a  simple 
form.  Instead,  we  will  express  these  two-layer  coefficients  in  terms  of  the 
half-space  reflection  and  transmission  coefficients  using  a  matrix  method. 

Let  the  amplitude  vectors  of  the  incident  and  reflected  waves  in  region 
0  be  a  and  b,  respectively,  as  shown  in  Fig.  3.  A  and  B  are  the  ampli¬ 
tude  vectors  of  the  upward  and  downward  propagating  waves  in  region  1,  respec¬ 
tively.  C  is  the  amplitude  vector  of  the  transmitted  wave  in  region  2.  Then 


these  amplitude  vectors  satisfy  the  following  matrix  equations 


**01 

lT01 


10 

*10J 


A 


(34a) 


B 

X 

fR 

K1 2 

c 

‘T12‘ 

(34b) 


where  Rg-j ,  ,  R-jg,  T^g,  R12,  are  the  half-space  reflection  and  trans¬ 
mission  matrices.  It  should  be  noticed  that  R^  and  also  account  for 

the  phase  difference  because  the  boundary  is  shifted.  Essentially  (34a)  and 
(34b)  reflect  the  boundary  conditions  at  z  »  0  and  z  =  -d,  respectively. 
From  (34b), 


B  -  R12A 


c  *  t12a. 


(35a) 

(35b) 


From  (33a), 


b  *  Rgia  +  T-jgB 


A  *  TqI a  +  RioB 


(36a) 

(36b) 


Substituting  (35a)  into  (36b),  we  express  A  in  terms  of  a: 


-*  -»  -•  ■> 


(37) 


A  *  (I  -  Toia* 

Combining  (35a),  (36a),  and  (37),  b  Is  given  by 

b  3  (R-jq  +  tioR12^  "  R10R12^  T01>  a* 

If  we  let 

b  =  Ra  (39a) 

A  =  Da  (39b) 

B  =  Ua  (39c) 

C  5  Ta  (39d) 

then  R,  D,  U  and  T  are  considered  as  the  reflection  matrix,  downward 
propagation  matrix,  upward  propagation  matrix,  and  transmission  matrix,  re¬ 
spectively.  Using  the  equations  (35)-(38),  we  obtain 

R  *  **01  +  T10R12^1  “  R10R12^  1  T01  ^40a^ 

D  -  (I  -  R10R12)-’  T0,  (40b) 

U  -  R,2(I  -  R,0R12)-’  T0,  (40c) 

t-T12(I  -  R)0R12)-1  To,.  <404) 


■-&8&G&6X, 


Mcy.tifjifayxTCr^.P';^  w^anr.  t-.1 r :  '.»  .*.:  y.  t- -■.  -’.  'A  T:  •*.;••-.  --  --.  y.  ~-.v~-..  •-■.-*•. 
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T!2  3 


eik?zd  R  eikUd 
e  k12oo  e 

iklzd  iklzd 

.  R,2oee 


!k"  d  -1k2  d 

e  x12oHe 

„1k?zd  .  -’k2zd 

X12oV  e 


e-<d  R  ;kizd" 

e  K12eo  e 

-<*  >kud 

•  R12eee  . 

e-1klzd  e-ik2zJ 

e  H2eH  e 

-ik®dd  -ik2  d 

p  V  p 

e  *12eV 


o  e 
o  o 


°Q  e 
e  o 


^eH 


*eV 


(42e) 


(42f) 


The  equations  which  all  these  half-space  coefficients  govern  and  the  expressions 
for  them  are  shown  in  Appendix  A.  Substituting  (42)  into  (40),  we  finally  ob¬ 
tain  the  expressions  of  two-layer  coefficients  in  terms  of  half-space  coeffi¬ 
cients. 


R$a  -  *01 8a  ♦  Xeo<Ll°o  +  "2eo)Xoa  +  X8o<Ll°e  +  Ve^ea 


+  WL2°o  ♦  Mleo>Xoa  +  XBe<L2°e  *  "|%H 


ea 


A6o  "  X6oL1  +  X8eL2 


%e  *  Vft  *  XseHl 


b8y  ■  xSo<L10r  +  "2eY>  +  XBe(L2°v  +  Mlev> 

X6a  *  Wh'o*  +  VJ  +  W^oa  +  Ml‘ea> 


where 


wl> V«  ■  i'.~r  ■' : 


(43a) 

(43b) 

(43c) 

(43d) 

(43e) 


m 


1  -  s. 


1  -  s. 


S  *  o  R  +  o  R 
oo  o  oo  e  eo 


3  +  o  R  a 

oe  o  oe  e  ee 


S  *  e  R  +  e  R 
eo  o  oo  e  eo 


S  *  e  R  +  e  R 
ee  o  oe  e  ee 


D  »  (1  -  S  HI  -  S)  - 

'  oo' '  ee  oe  eo 


(44a) 


(44b) 


(44c) 


(44d) 


(44e) 


(44f) 


(44g) 


(44h) 


(441) 


6,  a  *  H  or  V 


r  ■  o  or  e. 


The  expressions  in  (43)  can  also  be  derived  independently  by  using  the  method 
of  series  expansion  where  we  express  each  two-layer  coefficient  as  an  infinite 
series  in  terms  of  half-space  coefficients  by  identifying  all  the  constituents 
and  summing  the  series. 


IV.  Reflection  and  Transmission  Coefficients 


First,  consider  the  wave  with  TE  polarization  incident  from  a  free  space 
(region  0)  upon  an  anisotropic  medium  (region  1)  as  shown  in  Fig.  4a.  The 
electric  fields  in  each  region  can  be  formulated  as  follows 


Eo<f>  *  V"koz)  e  +  R01HH  hVkoz)  e  +  'VlHV  vo(koz}  e 


z  >  0 


(45a) 


.  .  iK,°  •  r  .  iR.e  •  r 

"  *Ho  0(‘klz)  e  +  *:ie  e(klz>  e 


z  <  0. 


(45b) 


The  boundary  conditions  to  be  satisfied  at  z  *  0  are 


z  x  EQ(r)  =  z  x  E1  ( r) 


(46a) 


A  A 

z  x  V  x  E0(r)  *  z  x  7  x  ^(r). 


(46b) 


Substituting  (45)  into  (46),  we  obtain  four  algebraic  equations.  Solving  them 
for  four  unknowns,  we  obtain 


XHe'-^-TiMk02  +  kozk°z>  s1n* 


k 

e  p 


where 


cos2*  kp2(kl2koz  -  ko2klP  +  si"2*{kl2<koz  -  ki°Hkx2  +  kQZk°2) 


+  ky2|clz(kl2  ’  ko2)}  +  cos  *  sin  *  Mk°z  +  koz)(ke°  ’  klz} 


^kp2  +  klzkoz^’ 


el  '  elz 


1 —  (kx2  *  k 2  +  k*d2  -  »2p£))[kx2  +  k  2  +  k®d2 


1/2 

j2u(e1  +  elz)]  ■ 


fd - L__  2  -22.  {k  2(k.2k 

D  (k  +  k?  )  k  p  T  02 
e  oz  lz  p 

+  ky(ko2kl°z  *  kl2kzklz)  sin  lp} 


ko2k?zd)  COS  * 


gd  3  <kx2  +  (k  cos  ip  +  k°z  sin  ^)2 > 


"olHH  ■  +  7  7"  V  C0S  *  +  kyk?z  s1"  *>  -  7“  kx s1"  *  r 


—  k?  sin  <p  +  — — —  (k  2k?^  cos  ip  -  k?2k  sin  ip) 

gA  k  k  12  U.kk  p  lz  lz  y 

%  p  oz  d  p  oz 


Second,  for  the  incident  wave  with  TM  polarization  [Fig.  4b],  the  electric 
fields  satisfy 


iR  •  r  .  i&  *r  -  %  iK0*^ 

E0(r)  *  v0(-k0Z)  e  +  **01  w  vo(koz)  e  +  R01 VH  ho^koz^  e 

z  >  0 


E-j(r)  *  Xy0  ®(“k°z)  ‘r+XHe^z)e 


ed\  'V  *  'r 


z  <  0. 


Applying  the  same  boundary  conditions  as  (46),  we  obtain 


U,  2k 


<kozkP2  cos  *  +  kykozklz  s1n  *> 

ue  p 

gH  2k  k  k.2k  k  -  k®J 
,  id  __o  Q2_i_ji  _oz - 15.  sin  <j> 


D«  k 
e  p 


k  +  ku 
uZ  KlZ 


"OIVH  *  7  7  V  cos  *  +  kyk°z  s1n  *>  ■  7"  kx  si"  *  Yf 

9d  p  p  d 

"om  "  -  yfrr  k°z sin  * +  777”  (k»2k'" cos  *  ■  k°^ sin  *] 

9d  p  oz  dp  oz 


d  p  oz 


Now  consider  the  wave  incident  from  an  anisotropic  medium  (region  1)  upon  an 
Isotropic  medium  (region  2).  For  the  ordinary  wave  incident  as  shown  in  Fig. 
5a,  the  electric  fields  satisfy 


E,(r) 


A  „  iRi  •  r  .  .  ik. 

°(-kl2)  e  +  Ri2oo  0(klz}  e 


i  R, 0  •  r 


r  ,  iK,e  •  r 

+  R12oe  •<kli>  e 


z  >  0 
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TT^TT-tt’ 


I 


E2(?)  ■  X120H  h2(-k2z>  e 


iR„  *  r 


o  .  „  iK,  •  r 

t  _t_  V  ..  t  U  \  ~  £ 


+  X12oV  v2(-k2z)  e 


z  <  0. 


Applying  the  same  boundary  conditions,  we  obtain 


9..  K  kT,  "  k. 


=  iy._i  Hz _ Hi 

12oo  a  c  k°  +  k 

9d  he  Klz  k2z 


where 


Ge  -  cos2*  kp2(k12k2z  +  k22k®“)  +  sin2i)>{k12(k2z  +  k*“)(kx2  -  k2zk°z) 


-  ky2klz(kl2  ‘  k22)}  +  cos  *  sin  *  ky(k2z  “  klz)(klz  "  klz)(kp2  "  k2zklz} 


9U  *  {kX2  +  ^ky  C0S  *  •  klz  sin 


1/2 


u  1 

R,2oe  *  —  —  2kxk“2(k2z  -  k°2){k  sin  *  +  k2z  cos  t)  sin  * 


9d  F 


where 


Fe  -  cos2*  kp2(k12k2z  +  k22klz^  +  sin2^{ki2^k2z  +  kIzHkx2  +  k2zklz^ 


eu. 


+  kv2klz(kl2  '  k22)}  "  cos  *  sin  *  ky(klz  +  k2z)(klz  +  klzNkp2  +  k1zk2z) 


U  *  i 
u  1 


£.  I  1 / 2 

(kx2  +  ky2  +  k^J2  “  ^pe^Ck*2  +  ky2  +  k-|z2  "  w2y(ei  +  elz^i 


i£l  ’  elz 


n-~.  V  . 


X1?«H  *  <kn2  cos  *  +  kvk?7  sin  ♦)  +  -^(k  2  cos  *  -  k  k?  sin  *) 
120H  9dkp  P  y  9ukp  p  ylZ 


k12kx  sin  * 


I  *  77-  -  Viz  s1n  *  -  -7^  kxklz  sin  *  *  TP  <kykll  sin  * 
k2zkp  l8d  %  Uu 

-  kp2kl“  cos  *)}  • 


For  the  extraordinary  incident  [Fig.  5b],  the  E  fields  satisfy 


-  .  ...  ilcT®  *  r  . 

^1 ( r)  *  e(kiz^  e  +  R12ee  e^klz^  e  +  R12eo  °^klz^  e 


.  V* 
0  \  .  1 


z  >  0 


^(r)  *  X12eH  h2(-k2z)  e 


iRn  •  r  -  iiC,  •  r 

!  ^  +  X12eV  vZ^~k2z^  e  ’  2  -  0< 


Applying  the  boundary  conditions,  we  get 


Uu  He 


UdFe 


where 


He  *  cos2*  kp2(k12k2z  +  k22k®z)  +  sin2*{k12(k2z  +  k®z)(kx2  +  k2zk°2) 


-ii 
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*  ky2klZ<kl2  *  k22»  -  cos  *  s1n  *  ky<k°2  +  k2z)(kl2  +  kU>V  +  k°2k22> 


*1260  ■  7T~~  kl2kx<kl2  ■  kl2)(k°z  *  k22)<ky  s1n  *  '  k22  cos  s1n  * 

ud  Fe 


,eu  ,.edw,,o 


12ev 


k„_k 


2z  p  tud 


—  (k?|k  sin  $  -  k  2kf^  cos  *)  -  k  k° 


‘lz^y  a'"  v  ~  *p_*lz  UU3  w/ - Viz  Sln  * 


*  ir  Vi* s1n  *  ■  kc2kiz cos  ♦> 


x  »  JL  j 

A12eH  k  ' 
p 


k  2  R  r 

-  —  k  sin  ^  +  -----  (k  2  cos  $  -  k  k°  sin  *)  -  -  12—  k,2k  sin  ^ 

,  U.  *  g  p  y  ,z  U  x 

{  d  su  u 


As  a  final  configuration,  the  wave  is  incident  from  an  anisotropic  medium  in 
the  lower  region  1  upon  an  isotropic  medium  in  the  upper  region  0.  For  the 
incident  ordinary  wave  [Fig.  6a],  the  E  fields  satisfy 


E,(r) 


*  n  iki  •  r  „  1 1\. 

°(kl2)  e  1  +  koo  °(-k?2)  e 


IR  °  •  f 


+  Roe  e<k?2)  e 


■d.  .V  •  ? 


z  <  0 


IK  •  r 

'oH  V"o2'  *  *  XoV  v0(k02)  e 


E0(?>  *  *„u  MO  ®  0  "  '  r 


z  >  0. 


Applying  the  boundary  conditions,  we  get 


3d  Ee  k°2  -  ko2 
3u  De  k°2  *  ko2 
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where 


Ee  -  cos2*  kp2(k,2koz  -  k02k^)  +  s1n2*{k,2(koz  -  k*d)(kx2  -  k^) 
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kozkp  9u  9d  Ud 


For  an  incident  extraordi nary  wave  [Fig.  6b],  the  l  fields  satisfy 
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Applying  the  boundary  conditions,  we  obtain 


_  ud  ie 

R  *  -  —  — 

ee  u  D 
u  e 


where 


Ie  *  cos**  kp*(k,*koz  -  k0*k*“)  ♦  Slr**{k1*(k0z  -  k*;)(kx*  -  k02k“z) 

+  ky2kl2(kl2  *  k02»  +  C0S  *  Sin  *  ky(k?2  +  k02)(kl2  '  k“><kp2  +  k02kl2> 


eH 


eV 


.h± 

U  1 

y 

u 

e 

,  JL. 

k 

' 

.  •  * 

p 

. 

k 

t  a  M 

0 

k 

k 

oz 

R 

,ed  ,.euw,.o 


1  R  R 

- k-.2k  sin  *  +  (k  2  cos  *  +  k  k?  sin  i|>)  — —  k,2k  sin 

li  1  *  „  p  y  iz  11  1  x 


9d 


—  (kizky  sin  ♦  -  kp2k*z  cos  *)  +  —  kxk°2  sin 
Uu  9d 


Ud  y 12 


c  2k?J 
P  IZ 


31 


V.  Conclusion 


The  dyadic  Green's  functions  (DGF)  for  an  anisotropic  (tilted  uniaxial) 
medium  have  been  evaluated  for  the  unbounded  and  the  two-layered  cases.  Inside 
the  uniaxial  medium,  there  are  two  characteristic  waves,  i.e.  an  o-wave  and 
an  e-wave.  A  wave  incident  with  either  TE  or  TM  polarization  produces  both 
TE  and  TM  waves  when  it  reflectes  at  the  interface  of  an  isotropic  and  a  uni¬ 
axial  ly  anisotropic  media  because  of  its  coupling  into  both  the  o-wave  and 
the  e-wave.  With  the  availability  of  the  DGF  for  anisotropic  media,  many 
layered-medium  problems  can  now  be  solved.  For  instance,  the  problem  of 
radiating  a  Hertzian  dipole  with  an  arbitrary  orientation  on  top  of  a  two 
layer  anisotropic  medium  can  readily  be  solved. 


r* 
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Figure  1  Geometry  of  an  anisotropic  medium. 


Figure  2  Geometry  of  a  two-layered  anisotropic  medium. 

Figure  3  Amplitude  vectors  of  the  incident,  reflected  and  transmitted 
waves. 

Figure  4  Reflection  and  transmission  from  an  isotropic  medium  upon  a 
anisotropic  medium. 

Figure  5  Reflection  and  transmission  from  an  anisotropic  medium  (upper) 
upon  an  Isotropic  medium  (lower). 

Figure  6  Reflection  and  transmission  from  an  anisotropic  medium  (lower) 
upon  an  isotropic  medium  (upper). 
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